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Abstract 

Coupled quantum harmonic oscillators, studied by many authors using many different techniques 
over the decades, are frequently used toy-models to study open quantum systems. In this 
manuscript, we explicitly study the simplest oscillator model - a pair of initially decoupled quantum 
harmonic oscillators interacting with a spring-like coupling, where the bath oscillator is initially 
in a thermal-like state. In particular, we treat the completely positive and trace preserving map 
on the system as a quantum channel, and study the truncation of the channel by truncating its 
Kraus set and its output dimension. We thereby derive the truncated transition amplitudes of the 
corresponding truncated channel. Finally, we give a computable approximation for these truncated 
transition amplitudes with explicit error bounds, and perform a case study of the oscillators in the 
off-resonant and weakly-coupled regime numerically. We demonstrate explicitly that the substantial 
leakage error can be mitigated via quantum error correction. 

PACS numbers: 03.67a, 03.65Yz 



I. INTRODUCTION 



One of the canonical physical models in quantum physics is that of quantum oscillators 
coupled with harmonic baths. The dynamics of such models and their variations has been 
extensively studied, using various techniques [THTT]. These techniques include Markovian 
master equations [12] . quantum stochastic processes and quantum Langevin equations [TlfTBT^ 
[T7] , Kossakowski-Lindblad equations [HI [19], methods in density-functional theory |20j, the 
standard techniques of perturbation theory, among many others [2TI [22] . 

Quantum channels [23J can be used to quantify the dynamics of a quantum system, 
and can be described by truncated transition amplitudes. In this paper, we approximate 
the truncated transition amplitudes of a given channel, where the truncation is performed 
with respect to the quantum channel's set of Kraus operators and its dimensions. A 
truncated quantum channel is a trace-decreasing quantum operation, and quantifies the 
partial dynamics acting on the system. Knowledge of the truncated quantum channel 
has utility - lower bounds on the performance of quantum error correction codes with its 
knowledge [2"5T[3"3"] . 

In this paper, we work towards quantifying the approximate dynamics of a pair of initially 
decoupled quantum harmonic oscillators interacting with a spring-like coupling, where the 
bath oscillator is initially in a thermal-like state. We work with a truncated subset of 
the model's Kraus operators, and thereby approximate its truncated transition amplitudes. 
We note that the Kraus operators of oscillator-bath models have also been approximated 
by various authors [HI EH [35] . Recently, Holevo also gave a formal exact expression for 
the Choi-Jamiolkowski operator for Gaussian channels [36], which describes the dynamics 
of coupled oscillators. Our contributions in this paper, are the explicit upper bounds on 
the approximation error of the truncated transition amplitudes of two quantum harmonic 
oscillators coupled via a spring-like interaction, where the approximation is an explicit 
summation of a finite number of computable terms, and dependent on the size of the input 
and output dimensions of the truncated channel. Our results can be used to explicitly study 
this toy model with rigorous error bounds. In particular, we numerically demonstrate and 
provide lower bounds for the leakage error, and show how this leakage error is mitigated via 
quantum error correction. 

The organization of the paper is as follows. In Section [TTJ we introduce the preliminary 
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material needed for this paper. In particular, we review L 2 (R) Hilbert space, quantum 
states, quantum channels, Hermite functions, and the linear canonical transformations for 



the quantum harmonic oscillator. In Section III, we give a treatment of the truncated 
dynamics of two quantum harmonic oscillators interacting with a spring-like coupling, and 
give explicit bounds on the error term induced by approximating the truncated transition 
amplitudes with a finite sum in Theorem III.1| In Section IV we give bounds on Hermite 



functions that are needed for the proof of Theorem |III.1[ Finally we apply our results 
explicitly in Section [V] in the case where the oscillators are off- resonant and weakly coupled. 



II. PRELIMINARIES 



In this section, we review the theory of L 2 (R) Hilbert spaces, quantum states and various 
representations of quantum channels, Hermite polynomials and functions, and coupled 
harmonics oscillators. 



A. The L 2 (R) Hilbert spaces 

We refer the reader to [37] and [38] for an introduction to separable Hilbert spaces. In this 
paper, all Hilbert spaces are complex and separable, have countable bases, and are typically 
isomorphic to the set of square integrable functions L 2 (R) which necessarily have infinite 
dimensions. Let W be the dual space of Hilbert space "H. For Hilbert spaces % and /C, let 
L(7-L,K,) denote the set of linear operators mapping 7-L to /C, and let L(7i) := L{T-L,T-L). Let 
03("H,/C) denote the set of bounded operators in L("H,/C). 

We use the Dirac's 'kef \ipn) to denote a function if) in the Hilbert space "H. For Hilbert 
space H, let {\jn)}jehi denote its generic orthonormal basis. We also use \if)%, (fie) and 
\if)n)\<PK) to denote \ifjy) £S> \<Pk). We denote the inner product of Hilbert space H using the 
Dirac notation as (-|-)% which is a sesquilinear form. For functions if>,(p G 7i — L 2 (R), their 
inner product in the position basis is (ip\ip) H := f R ip{x)*ip(x)dx. The physicist's position 
and momentum operators x% and are sesquilinear forms that map the tuple (\if>n)i I'Ph)) 
to f R ip(x)*x(p(x)dx = (4>h\xh\¥h) and f R ip(x)* ^(p(x)dx = (iPhIphIvh) respectively. We 
also denote the identity map on % as 1%. 
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B. Quantum states and channels 



We refer the reader to [23J for an introduction to quantum states and channels. Define 
the set of quantum states on Hilbert space % to be ®(H), the set of all positive semi-definite 
and trace one operators in 53(H). When p G 3)(H <S> JC), we denote the partial trace of p on 
Hilbert space H as Tr w (p) := (jn\p\jn)- 

A quantum channel $ : 53(H) — > 53 (/C) is a completely positive and trace-preserving 
(CPT) linear map, and its non-unique Kraus representation is [39W4"T] 

$(p) = ^K P Kt, J2 KiK = 1 n 

Ke.s Ke£ 

where A C 53(/C,H) is called the Kraus set of $. We denote the basis-dependent matrix 
elements of the Kraus operators by Kjj/ so that for all Kg^, 

3,3' 

We define the transition amplitudes of <3? with respect to the Kraus set A to be 

rr^'^E^'^ (!) 

a sum of the product of two Kraus operators over the entire Kraus set. Now let p G 2)(H) 
and $(p) G 2)(/C) have the decompositions 

P = ^2pa,b\b H )M, $(p) = $^Pa',&'l^)«;l 

a, 6 a, 6 

so that in the Kraus representation, 

(b' K m P )\a' K ) ^ (b' K \j2^\o 

KeJi 

= ( b 'ic\ Yl ^ K j,A?K)(jn\^PaA h n)M ^2 K l,k'\ k n)( k '/c\ a '/c) 

Ke£ a,b k,k' 

= YYYY K 3d' K kMP^ b '\j')ic(j\b)n(a\k)n(k'\a')ic 

Ke£ a,b k,k' 

= I Kfe ' 6,K a,a' ) P a , b = r 4>?i ^ ^ P*fi- ( 2 ) 

a,6eIM \KeJ? / a.belM 

Hence 7^ a ' 6 )~*( a ' b ) quantifies the transition amplitudes of (6«|p|a«) to (6^|$(p)|a^). For the 
purpose of quantum information processing, it may not be necessary to work with the full 



Kraus set .ft. In this paper, we instead restrict our attention to the truncated Kraus set fi, 
which is some appropriately chosen subset of the full Kraus set. This truncation procedure 
approximates the channel well if the truncated Kraus set comprises of the 'typical' Kraus 
effects. For the purpose of quantum error correction, partial knowledge of the channel 
is already of great utility, and recovery channels can be constructed based on this partial 
information to give lower bounds on the entanglement fidelity of specifically chosen quantum 
codes. Hence in this paper, the truncated transition amplitudes j^ ' 6 ) - ^ ' b ) pj a y a central 
role in quantifying the truncated dynamics of the channel 



C. Hermite polynomials and functions 



For n e N, define the Hermite polynomials H n (x) and the Hermite functions ip n (x) to be 



2 d n 



H n (x) := (-l)V — e 



ip n (x) :-- 



e l x *HJx) 



(3) 



For example, Hq(x) = 1 and H\(x) = 2x. The reader can refer to [12] for the properties of 
the Hermite polynomials and functions. For c > 0, also define the rescaled Hermite function 
to be 



1pn,c{ X ) = ( x H\lpn,c,H) ■= \fc^ n {cx). 



(4) 



As stated by Watson 
that is 



s ^z n ^ n (x)^j n {y) 



Mehler's formula applies in the case when \z\ < 1 and z is real, 

Axyz — (x 2 + y 2 )(l + z 2 ) 



n=0 



Trfl - z 5 



cxp 



2(1 -z< 



(5) 



Mehler's formula also holds for all complex numbers \z\ < 1, with the series converging 
uniformly and absolutely [Theorem 23.1 



D. A pair of harmonic oscillators and their linear canonical transformations 



1. The classical model 



We refer the reader to jl5] for an introduction to the quantum harmonic oscillator. Define 
the classical Hamiltonian of a classical harmonic oscillator with mass m, resonant frequency 
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u), position coordinate x and momentum coordinate p to be 

Hmuysp ■= ^— + \mu 2 x 2 . (6) 
2m 2 

The model we study has the classical Hamiltonian 

where Hi Qtt0 := — y) 2 is the classical Hamiltonian representing the spring-like interaction 
between the oscillators where > 0. The spring-like interaction H int o introduces quadratic 
terms ^ and ^- into H, effectively renormalizing the oscillator frequencies from w x o and 



w yj0 to u x := Y^io + ^ and w y := y w y,o + ^ respectively Hence when F int := -fory, 

In an experimental setup, it may be impossible to turn off the interaction between the 
two oscillators. Then the physically measured oscillator frequencies correspond to the 
renormalized frequencies. Therefore, we work with the renormalized representation of the 
model Hamiltonian given by ([7]). 



2. The quantized model 

Define the Hamiltonian of a quantum harmonic oscillator with associated Hilbert space 
H, mass M > 0, resonant frequency u > 0, position operator x% and momentum operator 
Pu to be 

Hm^,^ := || + \Mu?x 2 H . (8) 

The set of rescaled Hermite functions {1^ /thj n,)}nm is the set of energy eigenfunctions 
of the Hamiltonian tt(M,u);x n ,p n )- Let the Hilbert space of the first and second oscillators 
be X and y respectively, both isomorphic to L 2 (R). Define x := <g> ly, y := lx <8> £y, 
Px '■= Px 8> ty and p y := 1^ ®py. Then the quantized model Hamiltonian (JTI) is 

H = H^^.^^.^ g) ly + 1^ £g> H myft;ijIi ^ + H int 

HrrijtjUjtjx^ai ^"lyi^y^iPa ^int (9) 

where H int := — kxy is the quantized interaction. 



The coupled quantum harmonic oscillators can be decoupled by a linear canonical 
transformation of the oscillator positions and momenta (H)]. Define the rotation matrix, 
the rotation angle, and the rescaled mass by 



R 



cos 8 sin 9 



sin V cos ( 



9 := - tan 
2 



-i 



2k/m 
uj 2 — u 2 

y x 



, m := y/m x m y 



(10) 



respectively. The use of straightforward trigonometry then gives 



cos 9 




1/2 



1/2 



1 



i , 4fc 2 /m 2 

1 + K^IF 



sin = -= I 1 

y/2 



1 



i I 4fc 2 /m 2 

1 + K-^) 2 



mi 



where 



4k 2 /m 2 
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fc 



W y,o - W x ,„) + 



m Y — m 



Note that the rotation angle quantifies the strength of the coupling, in the sense that cos 9 ~ 
1 and sin# w when the coupling constant is small and the oscillators are off-resonant. 
Define the normalization parameter /j, := Then we choose the transformed position 

and momenta operators to be given by 



(u,v) : = R(/x 1 x,fiy), 



(12) 



where (xi,x 2 , ...) denotes a column vector. The quantized Hamiltonian is H = H m ^ u;ii ^ + 



H m , u „ ; i,p» where 



\ 





II 


' cos 2 01 


1 -- 






l sin 2 # j 


m 



\ 




+ — sin(20). (13) 
m 



Note that the frequencies o; u and u v are real as long as the original oscillator frequencies 
u Xj0 and uj 7)0 before renormalization are also real, because the renormalized frequencies u x 
and ujy increase as the coupling strength k increases. 

The linear canonical transformation that decouples the pair of harmonic oscillators is 
not unique. We chose the transformation that gives the same mass m for the decoupled 
oscillators, so that the only parameter different between them are the frequencies u u and 
u v . Since the transformation we have performed is canonical, [u,v] = \p u ,Pv] — [u,p v ] — 
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[v,p u ] = 0, [u,p u ] = ifo and [v,p v ] = ih. Hence there exist Hilbert spaces U, V isomorphic to 
L 2 (R) such that X ®y = U ®V , u = x u ® ly, v = t u <8>£y, p u = Pu® H-v an d — % ®Pv- 
Let c u := and c v := ^^p, c x := y^f^ and c y := Then the set of 

eigenstates of the uncoupled Hamiltonian H mx ^;x,p x + H my;U)y; y^ y and the full Hamiltonian 
H are {|^ K , Cu ,w, ^ X)Cv ,v)} k , x gn and {|^,c x ,*, ^/ 1 c y ,y)}j^eN respectively. 

III. TRUNCATED DYNAMICS OF THE INTERACTING SYSTEM 

This section highlights the main results of our paper. We provide a computable 
approximation to our physical model's truncated channel with corresponding error bounds 
that are simple to describe. 

A. The general model 

The Hilbert space of our model has the general form T-L = X <g> y where X and y 
are separable Hilbert spaces of the system and the environment respectively. Our model's 
Hamiltonian is 

H := H x ® t y + t x <g> H y + H int 

where H,H x ,H y and H int are (typically unbounded) Hermitian operators in the sets 
L(H),L(X),L(y) and L(H) respectively. The Hamiltonians H x and H y describes the bare 
dynamics on X and y respectively, and H int describes the system-bath interaction. 

Let the initial state of the entire model be p au := p ® 0y, where p £ ®(^) an d the bath 
state is a y = ^2g e ^Pe\^y) (^y\ £ ®(3^)- Let the time evolution operator of the entire model 
at time t be the unitary operator U t G 33 ("H). Then the time evolved state of system X at 
time t is 

pt := $ t (po) = Tr y U t p all U\ 

= £ (^l u *IJy><Jyl^®P^y><^l^<Jy|u!Ky> 

= £ (^|Ut|^Wo(^|Ui|^). (14) 
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Using Q, a feasible Kraus set and transition amplitudes for $ t are 



ft t :={^/pe(ty\V t \£y):£,£' eM} 



T, 



(a,b)^-(a',b') 



p e (b' x ,£ , y\U t \b x ,£y)(a x ,£y\Vl\a , x ,£ > y ) 



ft 



(15) 
(16) 



The transition amplitudes (16) of the full quantum channel may be impossible to evaluate 



because of the inifmite summation that is required. In view of this, we can instead evaluate 
the truncated transition amplitudes by truncating the infinite summation. These truncated 
transition amplitudes are the transition amplitudes of the truncated quantum channel. 



B. Coupled harmonic oscillators 

The approximate dynamics of coupled harmonic oscillators are still actively studied [13, 
[8] . In our paper, we use the model as described in Section II D 2 Let z u := e~ luJut and 



Z v . 6 



-liOyt 



For j G N, define |j x ) := \ip jjC ^x), \j y ) ■= \^j,c y ,y), \ju) ■= and 



:= iV'i.ov.v)- Then the unitary operator U f has the spectral decomposition 

U( ^ ^ ^/ Z U Z V Z U Z^ | K u , Xv) (^u? Xv | • 



(17) 



Let r = exp( 



k B T> 



G [0, 1), where ks is the Boltzmann constant and < T < oo is the 



effective temperature of the bath. The state of the bath with a Boltzmannian distribution 
is 



Thus pi = r (1 — r) in equations (15) and (16) 



(18) 



1. Kraus operators and transition amplitudes 



Using (15), the matrix elements of our Kraus operator K G At indexed by £, £' G N are 



(19) 
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The goal is now to find an expression for the truncated transition amplitudes for small values 
of a, b, a' and b' . We first give an expression for the transition amplitude with respect to the 
full Kraus set &t, which is 

= E ^(l-r)(^4|U t i6 x ^ y )(a x ,£ y |Uj|a' x ,4) 
e,t>ebl 

= r e (l-r)(b' x ,i' y \ ^ H ^ H K, Xv )(^X v KJ y ) 

x(a x ,£ y \ Zu K '~h; x '^\K' n , x' v )W u , x'yK, Q- (20) 

The matrix elements in the expression above can be simplified by expressing them in 
the x and y coordinates of the original oscillators. In particular, the expression above 
becomes an integral of the product of rescaled Hermite functions. To simplify notation, let 
u x , y '■= c u(-^rCOs9 + ^sinO) and v X)V := c v (— ^-sin^ + ^cos#) denote the coordinates 
of the decoupled oscillators in the basis of the original oscillators. By making appropriate 
substitutions, we have that 



u xy = \ x cos v + a / y sin 6 

w m x uj x y m y L)y 



m u v . m x u v 

v X y = —\l — — ismS + ,/ ycosu. (21) 

w m x uj x y m y u y 

The summation indices in the transition amplitude corresponding to the full Kraus set &t m 



(20) are I and £' respectively. In this paper, we choose our truncated Kraus set to be fi^t, 
where only the summation over £' is truncated. 

By applying Mehler's formula ^ on the variable £, the expression for the truncated 
transition amplitude is 



C™=S £ C*^/M,«^*<X.Xl^Jm. (22) 
£'<l K,xeN x y V v / 

where the path-dependent and time- independent transition amplitudes 

f[a,b,a',b';£',K,K',x,x') are 

1 + r 2 ( 2 4ry 3 y 4 , 



f 1 + r 2 ( 

f[a,b,a ,b;£ ,k,k ,x,x] ■= / dxdyexp - — — . 

ix.veR* .. 2(l-r 2 )V 



:,yeR 4 

X (^a'^l)^' (^2)^6(^3)^(^4) ) (#(l/l)#(2/2, 



(23) 
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where x = (xi, X2, X3, X&) and y = (2/1,2/2, 2/3> 2/4) are the rescaled position coordinates of the 
system and environment oscillator respectively. The integral / can be more easily evaluated 
if we express it as a product of three integrals, in the sense that 



/[a, 6, a', £ , k, k , x, X ] — h' ,e' y , x ' h' /' ,K, x Jb,a,K, x ,K' , 



where the integrals are 



Ia',t',K',x' - z 



R2 



@x,s,(a';f)V' K '( M x i2/ )V ; x'(^^) da; dy 



(24) 



(25) 



b,a,K,x,K 



®w,x,y,z 

(b, a; r)i! K (u Wi y)i! x (v W) y)i!^(u XyZ )ifj x/ (v XyZ ) dw dx dy dz., (26) 



and the kernels are 



^w,x, y ,z{i,3'^) '■= ipi{w)ipj(x) exp 



1 + 



2/ 



4r 



r 2/z + ^ 



(27) 
(28) 



2(1 -r 2 ) V l + r 2 ' 

Now our truncated transition amplitude is still a sum over an infinite number of integrals, 
and hence we intend to approximate it by the finite sum 



Aa,b)^{a',b') 
L,N,t 



E zr K 'zr x 'f[a,b,a',b';l' 



1 



,K,K,X,X 



i'<L k,x<N 
k',x'<N 



Lo x oo y y 7r(l + r) 



(29) 



for some positive integer N. In Theorem III. 1 , we prove that the absolute value of the error 



term by approximating the truncated transition amplitude (22) with (29) vanishes as N 



becomes large while a,b,a',b' and L remain small. The proof of our theorem uses mainly 
the Cauchy-Schwarz inequality, Lemma IV. 3 and Lemma |IV.4 The intuition behind our 
technical lemmas are elementary consequences of the behavior of order n Hermite functions 
in the oscillatory interval [—\/n, y/n\ and outside of it. Within the oscillatory interval, 
Hermite functions have amplitudes that vanish as n gets large. Outside of the oscillatory 
region, Hermite functions have exponentially small amplitudes as their arguments becomes 
large. Thus we construct our upper bounds for the integral of the product of Hermite 
functions by performing the integration separately in two overlapping regions, as depicted 
in Figure [TJ 

Theorem III.l. Let m x , m y , w x , u y , t > 0, k > and < r < 1 be real numbers, and D,L 
and N be positive integers. Let C be a constant that depends on m x ,m y ,co x ,Uy and k (see 
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(34)), and A,A,B,B be constants that depend on D and L (see (35)). Let T { 



-,(a,b)->(a' ,b') 



be 



the truncated transition amplitude defined in (22) have approximation A^'^^ a ' b ' given by 



(29). Then for all integers < a, 5, a', b' < D, 



(a,b)^(a' ,b') 



n, 



A 



(a,6)->(a',6') 



L,N,t 



aVaab 



-JVC/2 



+ 



-JVC 



+ 



\Ja^b 



-JVC 



9 (AT _ 1)3 3JV 5 /4(AT - 1)3/2 I _ e -C/2 ^5/2 ^ _ e -C/2)2 1 N 5/2 ^ _ e -C/2)2 

2 



N 5/4 ( 1 _e-c/2)( 1 _ e -C) 1 '""( 1 _ e -c) 2 y ^ tt(1 + r) ' 

Remark III. 2. Observe that when C is a very large number, the upper bound of the above 
theorem is dominated by the expression 

AA 2 B 



< 



14.7103(L + l)n%n 2 L / uj u uj v 1 



(AT -1)6 



CO x U! y 



1 



81(jV — I) 6 \ u x UJ y ]j 7r(l + r) 
as jV becomes large (and hi is defined as (max <j<i H^Hi)- 

Remark III. 3. If the parity of a + b differs from that of a' + b', the approximate truncated 
amplitude is necessarily identically zero for all positive integers L and N. This is a result 



of a simple parity counting argument after noting that the I-type integrals (25) and J- type 



integrals (26) are zero whenever the parity of the sum of their indices are odd. Hence parity 



is conserved with regards to our physical model. 



Remark III. 4. The bounds of Theorem III. 1 can be substantially tightened using information 



pertaining to the I-type integrals (25), which are substantially simplier to evaluate than the 



J-type integrals (26). Using bounds for the J- type integrals (37), we have that for integer 



N' greater than N, 



T (a,b)^(a' ,&') A (a,b)->(a',b')\ ^ V-^ . r j |d/„J vv 



,(a,&)->(a',b') 



A-5/2 



0<e<L N<n,n'<N' 
N< X ,X'<N' 



+ 



\T, 



(a,b)->-(a' ,b') 



A 



(o,6)->-(a',6') 
L.N't 



(31) 



Proof of Theorem III. 1 Our goal is to obtain upper bounds on each of the integrals / and 
J defined in (25) and (26). Applying the Cauchy-Schwarz inequality on \I a' ,1' ,k' ,x'\ §iv<3s 



\Ia>,i\K',x'\ — 



Q x ,y{a';£')\ip K i(u x , y ) 2 dx dyA \<d Xyy (a';£')\ijj xl (v X) y) 2 dx dy. (32) 

R 2 ' V JR 2 
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We similarly use the Cauchy-Schwarz inequality to obtain an upper bound of the absolute 



value of ( 26 ) , which is 



\Jb,a,K,x,K', X ',r\ <\l J \Qw,x,y,z( b ^ a ^ r )\^^( U ^,y) 2 ^K'(u XtZ ) 2 dx dy 



x 1/ / \Q w .x,y.z{ b ^ a ^ r )\i J x( v w,y) 2 ' l Px'( v x,z) 2 dx dy. 
Ir 2 



For the purpose of using Lemma IV. 3 and Lemma IV. 4, define the constants 



Ci = mm 



— - — cos 9, 

m x u x 



sin# 



m y u y 



■'" u cos 9 



c 2 = mm < , I — - — - sin 9, 

m x u x 

C = nun{l/(4^),l/(4^)} 



m y uj y 



(33) 



(34) 



and 



A = 1.74 2 max 



(max ||^||i) (max ||^||^ , A = {A.lA){2e^) D+L ^ D\D^> L\L^ 



B = 57.6 ( max 

0<j<D 



, B = (39.Q)(2e- 1 ) D D\D D 



(35) 



Noting that sin# and cos# are positive by definition (see (11)), and using the definitions of 



Ci, C2 and C with Lemma IV. 3 we have that 



I a >,e>, K ',x>\ < yJA(K')- 5 / 2 + Ae- ,cl JA(x')- 5 / 2 + Ae-x' 



>c 



(36) 



and we have a similar upper bound of \Iv,z\k,x\- Using Lemma IV.4, we have that 

\\a, w >,x'A < Jb(kk')-*/ 2 + Be-(K+«')cJB(xx')- 5/2 + Be~(x+x')c_ (37) 



By expanding out the terms of the products of the upper bounds given by (36) and (37) 



our upper bound on the absolute value of ( 24 ) 



f[a,b,a',b'-l',K,K',x,x'] <^W^W- 



where 



A 2 B 



k 5 (k') 5 k 5 (k') 5 / 2 
A 2 B 

+ K 5/2( ft ./)5/2 e 



AAB , c AAB kC A 2 B 

e + rnr, H 



-(k+k')C 



-f/.+«')c , AAB c _( k+2k >)c 1 AAB c _ (2k+k / )c ^ £2g e -(2K+2K')c 



,5/2 
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By the subadditivity of the square root function, we have that 

V - K 5/2^/)5/2 ^ K 5/2^/)5/4 C ^ ,.5/4(^)5/2 C ^ K 5 /2(«')5/4 

( K+K ' )c/ 2 "SAAB -( K +2Kf)c/2 , }/aX§ , 2K+Kf)c/2 I 7 2 5 

+ K 5/4 (K0 5/4 e + K 5/4 e + (ft;0 5/4 e 

The summation of the above expression over re and re' can be seen as an inner product of 
vectors with exponentially decaying terms and polynomially decaying terms respectively. 
We hence apply Holder's inequality for sequence spaces on the summation of the above 
expression over k and k'. We thereby obtain an upper bound of the sum in terms of the one- 
norm of the vector with exponentially decaying terms, and the infinity-norm of the vector 
with polynomially decaying terms: 

2^y W ^'^ VAB {l i X ) + AT5/4 {2^ X ) 1 _ e -C/2+ N5 /2 (l_ e -C/2)2 



>7V X>/V X>7V 



e -^ c 2\/a~AB e~ 3NC / 2 e~ 2NC 

+ , TE /o 77 77777777 + 777T77 77 77777777 777 + V A 2 B - 



N 5/2 (i_ e -c/2)2 /V 5 / 4 (1 -e- c / 2 )(l -e- c ) (l-e^) 2 ' 

Now the integral x~ 5 ^ 2 dx is an upper bound of the sum X}x>tv A~ 5//2 by the convexity 
of the integrand. Hence ^x>/v^~ 5 ^ 2 — f(^ — |)~ 3 ^ 2 > an<1 we can subsitute this bound into 
our upper bound of the square of J2 K k'>n summed over the index < £' < L to get 
the result. □ 



IV. BOUNDS ON HERMITE FUNCTIONS 



This section provides the main technical lemmas that are used to obtain error bounds on 
our approximation to our truncated transition amplitudes. The main technical tools that we 
use in this section are Alzer's sharp bounds on the gamma function [49J and bounds on the 
Dominici's asymptotic approximation of Hermite functions with error estimates by Kerman, 
Huang and Brannan [50] . 

Lemma IV. 1. For all positive integers n and reals x G [— y/n, y/n\, we have \ip n (x)\ < 
1.74n~ 5 / 4 . 

Proof. This proof combines Alzer's sharp bounds on the gamma function |39] with uniform 
bounds on the envelope of the Hermite functions in the oscillatory region by Kerman, Huang 
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and Brannan [50]. Using Kerman, Huang and Brannan's result (see equation (2.1) and (1.4) 
in [SO]), for x G [—y/ri, ^Jn\ we have 



<2- 3/4 v / 35 



-n . 



V2T((n/2) + 1) 

Alzer's sharp bounds for the gamma function are that for all n > 0, 

r(n + l) 1 1 5 

1 < ^ < H rf. 

v^^Hnsinhi)" 72 1620 

Note that for real n > 1, we have 1 < (nsinh ^)™^ 2 < 1.085. Hence we have that for positive 
integers n, 



Jr( n + f) A /V2^(^) n (1.085)(l + jj^) 

^0 < V ^i±_A ™> < 0. 9 30 8(2 »/>-V, 

Hence Kerman, Huang and Brannan's upper bound on the envelope of the Hermite function 
for x G [—y/n, y/n\ and n > 1 becomes 

|V>„(*)| < 2- 3 / 4 v^ ^ 1/ ^"" /2 ^ + + 1) 1) n- 1 < 1.74^. 

□ 

The next lemma provides a rather coarse upper bound on the absolute value of the 
Hermite function, with maximum utility in the monotonic region of the Hermite function. 

Lemma IV. 2. For all reals \x\ > 1 and integers n > 0, we have \tp n (x)\ < 2 n yj ' ^r^=e~ x2 / 4 . 

Proof. Using the Maclaurin decomposition of the Hermite polynomial H n (x) jSTJ, we get 



e -x 2 /2 0"f)l o, 

\M X )\ < n\n{2\x\) n ^== = x —^n\x\ n e~ x /2 . 



It is easy to verify that sup{|x| n e x ^} = ( — ) . Hence when \x\ > 1, 



2n 



n/2 



\^n{x)\<2 r 



□ 
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Lemma IV. 3 provides upper bounds on the one-norm of the product of Hermite functions 
in terms of the order of the Hermite functions, and is crucial in obtaining upper bounds on 
our error estimates. 



Lemma IV. 3. For real numbers a,& ^ 0, let c = min(|a|, |6|). Let j,k,n G N and n > 1. 
Then 



R' 2 



Max + byf\^(x)My)\ dxdy < O-^hmn umni 



+ (4.7A)2 j+k \r — 



e j+k 



Proof. We split the region over which the integral is performed into two overlapping regions 
A\ and A 2 (see Figure [T]), where 



A\ = {(x, y) : \ax\ + \by\ < \/n} 



and 



n 



A 2 = {(x,y) :x 2 + y 2 > ^=-} 

V 2c 



Then using the Charlier-Cramer bound [51] which states that sup na . 1^(^)1 < J// , and 




FIG. 1: The area inside the square corresponds to the region Ai (where \ip n (ax + by)\ < 1.74n -5 / 4 ), 
and the area outside the circle corresponds to the region A 2 (where all Hermite functions decay 
exponentially) . 



the uniform upper bound of the envelope of the Hermite polynomial in the oscillatory region 
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as stated in Lemma IV. 1 we get 



ip n {ax + by) 2 \ip j (x)ij; k (y)\ dxdy 



< 



1.74 2 



The integral of \ij)j{x)il)k(y)\ over the region Ai is at most 

(1.74 2 )||^||i||^l|i 



5/2 



\i*j(x)i>k(y)\ dxdy + 



1.086435 2 



7T 



\ipj{x)ipk{y)\ dxdy. 



n 



5/2 



Using the exponential upper bound of Lemma |IV.2[ the integral over the annulus region Ai 
is at most 

1.086435 2 2i+fcj jj\k\jik k 



< 



'it 

1.086435 2 



2> 



i+k 



71 



j\k\jik k 



-(x 2 +y 2 )/4 



dxdy 



A 2 



2tt 



re 



- 2 / 4 dr 



r> 



\ 2r 



.086435 2 2^ fc+ y^^(2e-"/( 8c2 )). 



Combining the upper bounds for region A\ and A2 then gives the result. 



□ 



The following lemma is needed to obtain upper bounds on the absolute value of 
\Jb,a,K,x,K',x'Ai an d 15 similar to the preceding lemma. 

Lemma IV. 4. Let 7^ a, b e R, c = min(|a|, |6|) and j, k, n, n' G N where n, n' > 1. Then 



tp n (aw + by) 2 ip n >(ax + bz) 2 \Q W)X , yiZ (j, k; r)\ dxdy 



< 



57.6||^||i||^ll 



- + (39.6)2 j 



■ +k j\k\jik k ( w+ „/- )/ ( 8c 2 ) 



( nn /)5/2 1 V—/- y ei+ fc 

Proof. We split the region over which the integral is performed into two overlapping regions 



A± and A2 (just as in the proof of Lemma IV.3), where 



Ai = {(w, x, y, z) : \aw\ + \by\ < yn, \ax\ + \bz\ < Vn'} 



and 



A 2 = {(tu, x, y, z):w 2 + y 2 > I—, x l + z 2 > -=-}. 

y2c V2c 
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Then using the Charlier-Cramer bound, Lemma IV. 1 and Lemma IV. 2 we get 

ip n (aw + by) 2 ip n '(ax + bz) 2 \Q WiXjy>z (j, k;r)\ dw dx dy dz 



R 2 



1.086435 4 



f 1.74 4 f 

- / 7 — e ^,wU) k > r )\ dw dx dy dz+ I 
J M (W) 5/2 J A , 

Now observe that we can express the kernel Q wx ,y zkh ^ r ) as 



®w, x ,y,zU' k;r)\ dw dx dy dz. 



(->,..-.,./,.:.(./• k;r) = e W+ z ^ exp 



1 _ r 2 

and hence the absolute value of our kernel has an upper bound that factorizes, in the sense 
that 

|e W)W 0',*;r)| < e~^ 2+z ^\^(w)Mx)\- 

Hence the integral over the region A\ is at most 

1.74 4 (2 7 r)||^]| 1 ||^i| 1 
(rzn') 5 / 2 

The upper bound on the kernel also allows us to find that the integral over the region A 2 is 
at most 

1.086435 4 



< 



7T 

1.086435 4 



-2 J 



i+k 



j\k\jik k 



-K+- 2 )/4 e -|(^+, 2 ) dw dx dy dz 



A-2 



7T 



3/2 



e j+k 



-r 2 



"d 



ri>- 



< 



1.086435 4 (27r) 2 2i+fc jj\k\jik k 



71 



3/2 



oj+k 



V2c 

4e -(n+n')/(8c 2 ) 



T2> 



V2c 



'"d 



r 2 



<(39.6)2^y^^e-( n+n ')/( 8c2 ). 

Combining our upper bounds on the integrals in the regions A\ and A 2 thereby gives the 
result. □ 



V. A CASE STUDY: OFF-RESONANT WEAKLY COUPLED OSCILLATORS 

In this section, we perform a case study of the dynamics of our physical model when the 
oscillators are off-resonant and weakly coupled. In agreement with the standard results of 
perturbation theory, we find negligible amplitude damping in the truncated dynamics of our 
system. Moreover, we show that leakage error is the dominant error process, and can be 
mitigated via quantum error correction. 
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A. Parameters of the physical model and the truncated channel 



The amount of truncation in our truncated channel A/ is quantified by the parameter 
L = 2, and the order of our approximation to the truncated channel is quantified by the 



parameter N = 6 (see (29) for the definition of the truncated transition amplitudes). We 
restrict the analysis of our truncated channels to a 4-level system by setting the parameter 
D = 3. 

Table [I] shows the parameters used for our physical model. We use the SI units. The 
output parameters are numerically computed using floating point numbers with 1024 bits of 
precision, and are shown up to ten decimal places. For our choice of parameters, the positive 





Input parameters 




Output parameters 




Other parameters and bounds 


m x 


10- 6 (kg) 


w x 


1000499.8750624610 (Hz) 


Wo 


1 


< 1.8827925275534299 


m y 


2 x 1(T 6 (kg) 


w y 


10000024.9999687501 (Hz) 


ll^i 


1 


< 2.124503864054394 




10 6 (Hz) 


Wu 


1001004.5438332595(Hz) 




1 


< 2.285324224284738 




10 7 (Hz) 


U} v 


10000025.0001779494 (Hz) 


I|V>3 


1 


< 2.410237758997186 


k 


100 


Ul 


10000025.0001779494 








r 





U 2 


0.0000031958 












Vl 


-0.0000451621 












V2 


1.0000000000 









TABLE I: We tabulate the important parameters of our physical model. Here, u\,U2,vi and V2 



are defined implicitly in the equations u x>y = u\x + U2y and v x , y = v\x + V2U (see (21 )). Also note 
that ||VVi||l := J R \rpn(x)\dx. 



constants AAB, A 2 B, A 2 B, AAB, A 2 B < 10 288 , and are negligible, and hence Remark 
holds. 



III.2 



B. Approximate dynamics of the truncated channel 



We numerically evaluate approximate truncated transition amplitudes (29) corresponding 



to the transitions within the lowest energy levels of the system. The evaluation of each 
such approximate truncated transition amplitude A^q^ ^° ,b ^ is a sum of 7203 terms for our 
choice of L = 2 and N = Q. We obtain the corresponding error bounds of our approximation 
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from Remark HI. 4 with N' = 15. The error of each of our approximate truncated transition 
amplitude is at most 0.00084 - a negligible amount. We plot magnitudes of the non-negligible 
truncated transition amplitudes in Figure |2j 

.e terms of A 2 6 i 

1 and < a', b' < 3 that satisfy the relation 

a' — a V — b 



Numerically, we find that the non-negligible terms of A^rI^^ ,b "* have values of < a, b < 



e Z. 



(38) 



2 ' 2 

The above relation holds for two reasons. Firstly, the conservation of the parity of all 



truncated transition amplitudes as stated in Remark HI. 3 implies that the parity of a + b 
equals the parity of a' + b'. Secondly, the other negligible transitions are in agreement 
with the results of using perturbation theory on off-resonant and weakly coupled harmonic 
oscillators. Also note that the negligible damping from the first excited state to the ground 
state of our truncated channel suggests that the coupled oscillator model is inconsistent with 
the phenomenon of amplitude damping [23] even at zero temperature. 

The Choi-Jamiolkowski (CJ) operator of the finite input and output dimension channel 
$ with Kraus set & is the linear operator 



X* := l K M K l = EE ( E K ^' K M' ) \3K>3n)(kcM' 

Ke£ j,f k,k' / 

where the linear map |-)) : L(H, K.) — > K, (g) H is a stacking isomorphism such that 

iE a ^>o«i» : = E a ^>b«>- 



(39) 



hi 



hi 



We now construct a quantum operation A to approximate the truncated channel M . 



Abbreviating our approximate truncated transition amplitudes as A aha i, 
our approximation to the CJ operator of our truncated channel is 



A 



(o,6)-»-(a',fc') 
2,6,i ' 



/An 



X 



0000 ^0002 Aqiqi Aqi 03 

A)020 Ax>22 ^1012 ^0123 

^ioio A1012 Ann ^1113 

\Al030 ^1032 ^1131 ^1133/ 



(40) 



where the bases labeled by the rows and columns are |0, 0), |2, 0), |1, 1), |3, 1) and 
(0, 0|, (2, 0|, (1, 1|, (3, 1| respectively. The first and second entries of our bras and kets 
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0.5 



L 



Absolute Values of Approximate 
Truncated Transition Amplitudes A^]^ 




/UHU)| 
*2,6i 



,(1,1)^(1,3)1 



t0«) 



FIG. 2: Absolute values of approximations to our truncated amplitudes are depicted for our model 
with parameters given by Table Q. 

correspond to the output and input Hilbert spaces of the truncated map respectively. We 
plot the absolute value of some of these non-negligible matrix elements in Figure [2] 

Let the spectral decomposition of \ have the form % = Y^t=i W^i)Qn\ where |Aj)(Aj| are 
orthogonal projectors and the eigenvalues \ are in non-increasing order |52j. For A« > 0, let 
Ai be the image of the inverse map of the linear operator |-)) acting on \f\iW). Numerically 
diagonalizing the matrix x shows that it has only one dominant eigenvalue. Hence we use 
the quantum operation A(v) := AivA[ to approximate the truncated channel M where A\ = 
v / AT(A;o|0)(0| + /c 2 |2)(0| + /c 1 |l)(l| + A;3|3)(l|) and |Ai) = k \0, 0} +fc 2 |2, 0) + 1} + fa\3, 1). 
Observe that for i, j e {0, 1, 2, 3} we have 

^(|*>01) = Ax (h\i) + k l+2 \i + 2)) (A^OI + k* j+2 (j + 2|). (41) 



1. Leakage error 

Leakage error, a dominant process of our model, occurs when low energy states transition 
into higher energy states within the quantum system. We derive lower bounds on the 

21 



minimum amount of qubit leakage in our system, when the quantum channel is $>t with 



Kraus set &t given by (15). For qubit leakage to occur, it suffices to have the strict inequality 

oo 

Leakage($ t ,p) := £<z*|$t(p)|i*) > 0, (42) 

i=2 

for some density operator p supported on the span of \0x) and \lx)- The above expression 
quantifies the amount of leakage from the qubit state space of our system. 

We proceed to obtain a strictly positive lower bound on Leakage^, p). If $i = Af + Ai 
for some completely positive maps Af and Ai, then the complete positivity of Af and Ai 
implies that Leakage($ t , p) > Leakage (A/ - , p). Hence it suffices to obtain a lower bound for 
Leakage (Af, p). For our application, Af is our truncated channel. 

Assume that the initial state of the system isp= ^(|0^)(0^| + |1^)(1^|), the maximally 

mixed state in the qubit space. Then the leakage of the truncated channel Af is at 
least i( T (o.o)-(2,2) +T (i,i)-(3,3)y Thug we haye that Leakage(A/ -, p) > l(|4°;°)-( 2 ' 2 )| + 

e is given the upper bound in Theorem 



III.l 



In view of the 



data given in Figure [2j the amount of qubit leakage can actually be quite substantial. In 
particular, when t = 5 x 1CT 6 , the amount of leakage of $4 is at least 0.4 which is substantially 
larger than zero. 

The large amount of leakage from our qubit state highlights the importance of accounting 
for transitions of low energy states into excited states in oscillator systems, and also the 
problem of using the lowest two energy eigenstates as a basis to encode our qubit. 



2. Quantum error correction versus no quantum error correction 

We consider a universe with four harmonic oscillators X\,X 2 ,y\ and y 2 . We identify 
the oscillators X\ and X 2 as 2- type oscillators and the oscillators and 3^2 as y-type 
oscillators, with their parameters given by Table |T| Assume that there are only Xi~yi 
couplings and X 2 ~y 2 couplings in our universe. Suppose that a maximally entangled two- 
qubit state is initialized in the X\ and X 2 oscillators supported on their two lowest energy 
levels. We assume that the oscillators and y 2 are initialized in the ground state. We 
obtain a lower bound on the fidelity of the time-evolved states when instantaneous, identical 
and independent recovery operations are performed on oscillators X\ and X 2 . Denoting 
our truncated channel and recovery channel on a single system oscillator by Af and TZ 
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respectively, our lower bound on the output fidelity is 

fc = (Q\(n®K)o(M®W(\$)($\)\$), (43) 



where |$) = (10^,0^) + \^x 1 , 1*2))/ v 2. We now proceed to obtain a lower bound on 
the entanglement fidelity of our time-evolved maximally entangled state, with and without 
recovery. The maximally entangled state on two qubits written as a density matrix is 

mm = \ £ (44) 

i,je{o,i} 

By linearity, the action of our truncated channels on the maximally entangled state gives 

w®Aom($\ = l Yl a/-(i*xj|)®a/-(k><j|). (45) 

i,je{o,i} 

Let us denote the error of the truncated transition amplitude of A^iji to be so that 

Ue{o,i} 

= \ Yl ( A ynii + e ynii)(^^ (46) 
M'e{o,i} 

jiji£{0,l,2,3} 
j2J2e{0,l,2,3} 

If we perform no recovery operation, 1Z is just the identity map X, and we have 

= ($|( J2 W\ixi)<Jxi\)®W\i*W*\)) I*)/ 2 - (47) 

\(*J)6{0,1} / 



Dropping the labels on the Hilbert spaces of our bras and kets, we can use (46) to find that 



fx = k(o,o\ + (MDCAf® AOI*X*|(|o,o> + |i, l)) 



7 ^ ] iA+jixh e ijiiji)( A ijhj2 + e ijhj2) 1*1) ^2) (jl ? J 2 | 



ij6{0,l} 
il,il6{0,l,2,3} 
i 2 J 2 e{0,l,2,3} 

> \(Kooo + ^111 + <oi + Aiow) - 8 ( 2e + e ') (48) 

where e < 0.00084. Note that A 2 0101 + A 2 W10 can be negative. 

The Barnum-Knill recovery operator TZ BK [26] and the Tyson-Beny-Oreshkov quadratic 
recovery operator [301 EE] are near optimal recovery operators defined with respect to a 
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quantum operation A and a state p, and are equivalent when A has only one Kraus operator. 
For our application, we study the Barnum-Knill recovery, with p = (|0)(0| + |l)(l|)/2 and 
quantum operation A approximating the truncated channel M. Our choice of p shows that 
our quantum information is encoded in the trivial quantum code (no encoding). For our 
application, the Barnum-Knill recovery operator which is also a quantum operation is defined 
as 

TZ BK (v) = A{{A l A\y 1 l 2+ v(A x A\)- 1 l 2+ A 1 (49) 

where (Ai-Aj) 1 ^ 2+ is the square root of the psuedo-inverse of the operator (A1A.J). When 
we use the Barnum-Knill recovery operation TZ BK , the fidelity of recovery is 

/bk = (1>\(n BK ® TZ BK ) o (AT® j\0(|*><*|)|*>. > A? - |A2l + |A 4 3l + |A4 ' - \(2e + e 2 ). 

(50) 

We plot lower bounds of the fidelity with and without Barnum-Knill recovery in Figure |3j 
and demonstrate that a fidelity of more than 60% is still possible in spite of the leakage 
error and the use of truncated quantum channels in our analysis. 

VI. DISCUSSIONS 

The system we consider is described by a quantum harmonic oscillator coupled through 
a spring-like interaction to another initially decoupled harmonic oscillator. We provide 
approximations to the truncated transition amplitudes of such a system. The converging 
error bound of such approximations is our main result. Properties of the integrals of products 
of Hermite functions lie at the heart of the proof. It is also worth noting that any ensemble 
of harmonic oscillators with spring-like coupling can be analyzed similarly. 

We also show numerically that in agreement with intuition from perturbation theory, 
when the oscillators are off-resonant and weakly coupled, amplitude damping is a negligible 
physical process. We also use our truncated channel representation to show that qubit 
leakage can be a dominant physical process, and how Barnum-Knill recovery can help 
protect a maximally entangled state stored in two oscillators each coupled independently to 
distinct zero-temperature harmonic baths, in the paradigm of off-resonant and weak coupling 
between the system and a zero temperature bath. 
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FIG. 3: Lower bounds on the fidelity of an entangled state without recovery fx (dashed line) and 
with Barnum-Knill recovery /bk (solid line) are plotted with respect to time. 
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